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Abstract 

The restricted three-body problem posses the property that some classes 
of doubly asymptotic orbits are limits members of families of periodic or- 
bits, this phenomena has been known as the "Blue Sky Catastrophe" 
termination. A similar case occurs in the restricted four body problem 
for the coUinear equilibrium point named 1/2- We make an analytical and 
numerical study of the stable and unstable manifolds to verify that the 
hypothesis under which this phenomena occurs are satisfied. 

1 Introduction 

Few bodies problems have been studied for long time in celestial mechanics, 
either as simplified models of more complex planetary systems or as bench- 
mark models where new mathematical theories can be tested. The three-body 
problem has been source of inspiration and study in Celestial Mechanics since 
Newton and Euler, in particular the restricted three body problem (R3BP) has 
demonstrated to be a good model of several systems in our solar system such as 
the Sun- Jupiter-Asteroid system, and with less accuracy the Sun-Earth-Moon 
system. In analogy with the R3BP, in this paper we study a restricted problem 
of four bodies consisting of three primaries moving in circular orbits keeping 
an equilateral triangle configuration and a massless particle moving under the 
gravitational attraction of the primaries. In 1933 in the Copenhagen observa- 
tory, Elis Stromgren and co-workers performed a numerical investigation of the 
periodic orbits of the R3BP for the value fi = 1/2 of the mass parameter. In this 
investigation we can find a family of periodic orbits called class g. This family 
originates form periodic orbits which Poincare called premiere sorte, [TS], [50] . 
From this simple beginning, the family develops many variations and the natural 
end of this family could not be stated by the investigators in the Copenhagen 
observatory. Stromgren wrote in his conclusions that the termination of this 
family had to be an asymptotic periodic orbit spiralling into the equilibrium 
points Li and L5. Several subsequent numerical explorations were necessary 
to confirm this conjecture. However, an analytical proof of this conjecture was 
given by J. Henrard [1^ until 1972 and later by BufFoni [7]. In those papers, 
Henrard and BufFoni give sufficient conditions to prove the conjecture in the 
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more general framework of analytical Hamiltonian systems with two degrees of 
freedom. These conditions amount to a transversality condition on the stable 
and unstable manifolds. In papers like [9] |3] there exist numerical explorations 
of several families of periodic orbits of the restricted four-body problem for 
some values of the masses, in these works we can find that some families be- 
come asymptotic to a equilibrium point named L2 as the value of the Jacobi 
constant at the equilibrium point is reached. So as in the R3BP, this numerical 
evidence requires to perform an analytical investigation to verify that the con- 
ditions for which the "Blue sky catastrophe" termination ocurrs are satisfied. 
This paper is organized as follows: In section 2 we state the equations of motion 
of the restricted four-body problem and we give some basic properties of the 
problem. In section 3 we analyse the linealization of the equations of motion 
at the equilibrium point L2, in particular we will analyse the eigenvalues of 
the linear equations. In section 4 we give a brief introduction to the normal 
form theory and an application to the restricted four-body problem will be per- 
formed. In section 5 we analyse the normal form obtained and we will verify the 
conditions given in [24] to prove a topologically transverse intersection between 
the stable and unstable manifolds. 



2 Equations of Motion 

Consider three point masses, called ■primaries, moving in circular periodic orbits 
around their center of mass under their mutual Newtonian gravitational attrac- 
tion, forming an equilateral triangle configuration. A third massless particle 
moving in the same plane is acted upon the attraction of the primaries. The 
equations of motion of the massless particle referred to a synodic frame with 
the same origin, where the primaries remain fixed, are: 
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where fc^ is the gravitational constant, n is the mean motion, pf = {x — Xi)^ + 
{y ~ ViY is the distance of the massless particle to the primaries, Xi, fji are 
the vertices of equilateral triangle formed by the primaries, and (') denotes 
derivative with respect to time t* . We choose the orientation of the triangle 
of masses such that mi lies along the positive a;-axis and TO2, are located 
symmetrically with respect to the same axis, see figure [T] 

The equations of motion can be recast in dimensionless form as follows: Let 
L denote the length of triangle formed by the primaries, x — x/ y — y/L, 
Xi = Xi/ yi — iji/ L, for i — 1,2, 3; M ~ mi + 7712 + the total mass, and 
t = nt* . Then the equations ([T]) become 
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y + 2x-y = -2_^fJ-t 3 — (2) 
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where we have used Kepler's third law: k'^M — ti^L^, and the dot (') represents 
derivatives with respect to the dimensionless time t and rf = {x — Xi)'^ + {y — yi)'^. 

The system ^ will be defined if we know the vertices of triangle for each 
value of the masses. In this paper we suppose := /is = ^2 then /ii = 1 — 2/i, it's 
not hard to prove that the vertices of triangle are given as function of the mass 
parameter fihy xi = VSn, yi = 0, X2 = -^^ItLzM^ = -i, 3:3 = _^/l0^zM^ 
ys — h- The system B can be written succinctly as 



x-2y = (3) 
y + 2x = ily (4) 
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where 

n{x,y,fi) l{x^+y^) + y] 

i—i 

is the effective potential function. 
There are three limiting cases: 

1. If /i = 0, we obtain the rotating Kepler's problem, with mi = 1 at the 
origin of coordinates. 

2. If = 1/2, we obtain the circular restricted three body problem, with two 
equal masses TO2 — — 1/2. 

3. If /i = 1/3, we obtain the symmetric case with three masses equal to 1/3. 

It will be useful to write the system Q using complex notation. Let z — 
X + iy, then 

z + 2k = 2^ (5) 
oz 

with ^ 

n{z,z,fi) = -Izp + U{z,z,fj,) 

where the gravitational potential is 
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U{z,z,fi) = V ^' 
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Figure 1: The restricted four-body problem in a synodic system 
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and Ti = \z — Zi\, i = 1, 2, 3. are the distances to the primaries. System ([5| has 
the Jacobian first integral 

2n{z,z,fi) -\z\^ = C (6) 

If we define P — Px + ipy, the conjugate momenta of z, then system (|3| can 
be recast as a Hamiltonian system with Hamiltonian 

H = ^\P\^ + Im{zP)-U{z,z,n) 

= ^(pI+pI) + {yPx - xpy) - U{x,y,fi). (7) 

The relationship with the Jacobian integral is H = —C/2. The phase space of 
([t]) is defined as 

A = {{z,P) eCxC\z^ z„i^ 1,2,3}, 

with collisions occurring at z = z.^, i = I, 2, 3. 

In the restricted three-body problem there exist five equilibrium points for 
all values of the masses of the primaries but in this restricted four-body problem 
(R4BP) the number of equilibrium points depends on the particular values of 
the masses. Figure [2] shows the Hill's regions for a large value of the Jacobian 
constant C consisting of a large exterior region around the primaries and three 
components containing them. As the critical Jacobian constant decreases, the 
evolution of the Hill's region is shown in Figure [2j The smaller value of C is 
just above the critical value where the Hill's region is the whole plain minus the 
positions of the primaries. 

A complete discussion of the equilibrium points and bifurcations can be 
found in [TU] , [22] , [H] , [1] , [H] . For the particular value of the mass parameter 
presented in [9], the R4BP has 2 coUinear and 6 non-coUinear equilibrium points. 
We use the notation shown in Figure [3] for the eight critical points. 

3 Linearization around the collinear equilibrium 
point L2 

In the former section it was shown that the number of equilibrium points in 
the restricted four-body problem depend on the value of the masses of the 
primaries. In papers like [H] and [5S] we can see that the stability of the 
collinear equilibrium points also depends on the values of the masses, in this 
paper we will perform a thorough study of the stability of the equilibrium point 
L2 (see Figure [3]) and its consequences. 

If we write L2 = (^,0) for the coordinates of the equilibrium point then the 
Hamiltonian function centred at L2 reads 

= ^iVi+yl) + {x2yi - Xiy2) - xiC- C/(xi,X2,Ai) 

with 

3 

U[XI,X2,IJ) = — 
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Figure 2: Hill's regions for a large value of the Jacobian constant (top-left). 
Hill's regions for critical values of the Jacobian constant. The Hill's regions 
of the last row correspond to a slightly larger value than the critical one for 
illustrations purposes. 



Figure 3: The eight equihbrium points for ji = 0.019. 



where rf = (xi + ^ — ui)"^ + {x2 — ViY' , i = 1, 2, 3. and Ui, Vi denote the positions 
of the primaries. Expanding through second-order terms, we obtain 

H ^^{yj+ vl) + {x2yi - xiy2) ~ ^ {U^.x^xl + U^^x^xl + ) (8) 

There are no hnear terms because the expansion is performed at a equihbrium 
point and the constant term has been omitted because it does not contribute to 
the system of differential equations. The term 11x^^x2 is zero because the equihb- 
rium point is coUinear. If we ignore the higher order terms, the corresponding 
quadratic Hamihonian gives rise the foUowing Hamihonian matrix 

/ 1 1 \ 

-10 1 

a 1 
\ 6-10/ 

Where a — Ux-ixi and b = Ux2X2i the characteristic equation of this matrix is 

P(A) = A-* + (2 - a - b)X'^ + {ab + a + b+l) 

using the relation Qxixi = 1 + Uxixi, ^^2:2x2 = 1 + Ux2Xo we can write as in the 
R3BP: 

P(A) = A^ + (4 - nx,x, - ^X2X2)><^ + nx^X^^X2X2- (9) 

The analysis made in |19j states that at L2 there exist such that for fi < fXh 
the eigenvalues are iiwi, ±za;2, for fj, — ^i, is (iio;)^ with multiplicity two, and 
for fj. > fib the eigenvalues are off the imaginary axis. The value /ifc can be 
determined numerically as a zero of the discriminant of 

obtained form ^ if we make i] ~ X^. This value is approximately = 0.0027, 
it must be clear that uj is in function of a and b and therefore of fi. As we 



said in section 1, there exists numerical evidence that some famihes of periodic 
orbits of this problem undergo a "Blue Sky Catastrophe" (BSC) termination, 
therefore this suggest to verify the hypothesis given in [16] and [7] to give a 
rigorous demonstration of these phenomena. 



4 Normal form at the equilibrium point 

Perturbation theory is one of the few ways to study the behaviour of a real 
nonlinear system beyond its linear approximation, there exist a lot of papers 
dealing with this problem, for example in [32], [33] j [IZ] [131 IS] it can be 
found a good introduction to the perturbation theory. However, we show briefly 
Deprit's algorithm in order to calculate a normal form at the equilibrium point 

Let if (e, x) be a Hamiltonian, X{e, y) be a change of variables generated by 
a function W[e,x) and G{e,y) — H{e,X{e,y)). Suppose that H, W, G all have 
series expansions in the small parameter e. The algorithm of the method of Lie 
transforms is a recursive set of formulas that relate the terms in these various 
series expansions i,e., let 

H{e,x)^H",{x) + Y,[^]H°A^) (11) 

n— 1 ^ ' 

^ / \ 

G{e,y) = H"o{y) + Y.(-.)HS{x) (12) 

n=l ^ ■ ^ 

W{e,x)^Y.(^)^-+^^=') (13) 

the functions Hj, i — 1,2, .., j — 0, 1, 2... satisfy the recursive identities 

i 



fc=0 

For example, to compute the series G(e, y) through terms of order we need 
to solve the first homological equation 

H^^H'^ + {HS,W^} (14) 
which gives the term of order e, then we must compute 

Hi = + {ii?, W^} + {ijo, W2} 
and solve a second homological equation 

H^ = Hl + {H^,W^} (15) 
which gives the term of order in this way we obtain 

G{e,y) = H°,{y) + eHl{y) + ^^Hl{y) + .... (16) 

We can transform the Hamiltonian (|8| in the Hamiltonian (11) if we make 
the symplectic scaling x = e{xi, X2,yi,y2) & with multiplier e^^ in the 
Hamiltonian ^ and collect the homogeneous terms. 



7 



4.1 Linear normal form 



If we ignore the higher order terms in ([8]) we get the following quadratic Hamil- 
tonian 

H =^{yf + yl) + {x2yi - xiyz) - ^ {U^^^^^xl + U^r^x^xl) , (17) 

this Hamiltonian can be written as a quadratic form H — Sx, with S — 
— J4A, J4 denotes the standard four dimensional symplectic matrix and A is 
the matrix defined in the previous section. We seek a symplectic change of 
coordinates x — Pz such that the Hamiltonian (17) can be written as H = 
S*z where S* = —J^B and B = P^^AP is the normal form of A 



B = 
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—UJ 




























[ 


e 


iO 


/ 



(18) 



where e — ±1, this value depends on the basis. In papers like [H] and [S] there 
are detailed discussions on how to find the matrix P, in particular we follow the 
algorithm shown in [8| to achieve this. The symplectic basis required depends 
on the decomposition A — 'S + N where S is a real semisimple symplectic matrix 



and N is nilpotent matrix given by A'^ = ^ — S. 
and Ld as in the R3BP these matrices looks like 



In terms of the quantities a, b 



-(3c^2 





2a 



3a; 



1 

1 



a 



6-1) 



- 6 + a(3w^ - 2) 




26- 



a + 6(3a;2 + 6-2) 



3^2 + a - 3 



-(3a;2 + 26 + a- 1) 



3u;2 





+ 6-3 
2a + 6 - 1 




(19) 



and 



v 





+ 2a + 6 - 1 
-(a2-6 + a(cj2-2)) 




N 



-(a;2 + 26 + a-l) -(tu'^ + a-S) 



26 



6(a;2 + 6-2) 










26 ^ 



1 

2^ 
a - 1) 

-(cj^ + 6-3) 
-(cj2 + 2a + 6- 1) 
1 

(20) 

We must observe that the structures of S and N are equal to those of the R3BP 
in the non semisimple case, more precisely, if we substitute the numerical values 
of a and 6 we can see that iV 7^ but iV^ = and it follows that the nilpotent 
index equals to 1; therefore we can build the desired basis as it was done in the 
R3BP. First consider the symplectic product on R^, {x, y) = x'^ J^y, second, we 
need a initial vector to build the basis and to determine the value of e. As in 
the R3BP a natural choose is the vector zn = , = with ei the standard 

^/\{euNe,)\ 

i*, it is not hard to see that {ei,Nei) = N31 and therefore 



unitary vector in . 



(zo^Nzo) 



(ei,iVei) 

1^31 1 



31 



311 
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the value of N^i is approximately —1.82 and therefore e = — 1 as in the R3BP. 
The desired basis has the form : zi — zq + ^^{zo,'Szq)NY,zq, Z2 — ^Szi, 
Z3 = eNzi, Z4 = ^EA^zi. Therefore the desired matrix is P = col{zi, ^2, 2:3, Z4), 
this matrix has the form 



P = 



/ Pii pi4 ^ 

P22 P23 

P32 P33 

\ P41 P44 ) 



(21) 



of course the entries of this matrix depend of the entries of E and N (therefore 



of a, 6). The Hamiltonian ( 11 ) under the change of coordinates z = Px becomes 



i/(e,z) = i/°(2) + 



CX3 

E 



(22) 



Hq{z) is the linear Hamiltonian associated with the normal form B. 



Hamiltonian (111, recalling the expansion given in 
H'^{x) contain the terms of order n 



4.2 The higher order terms 

In this section we deal with the higher order terms in the expansion of the 

we see that the terms 
2 of the Taylor series of the potential 
U{xi,X2, 1^), therefore such terms have the form i?i(x) = a^xf + 63X1X2 + 
CsXiCcI + (^3X2, -ff2 (x) = 040:^ + 64x^X2 +C4xfa:| + ^4X1X2 + 64X2, ...etc. Here the 
coefficients of the these homogeneous polynomials correspond to the respective 
derivatives of order n + 2 in the expansion of U (xi , X2, /i) evaluated at the origin. 
Because the equilibrium point is collinear we have that ^^{xi,0, fj.) = and 
therefore g^mg^^ U{xi, 0, /i) = for all n. A straightforward calculation shows 



dHJ_ 
dxl 



= M 



9(x2 - 1/2) 15(x2 - 1/2)= 



9(x2 + 1/2) 15(x2 + 1/2)= 



-(1 - 2^l)^, 



9X2 



15xi 



if we evaluate in X2 = we see that r2 — ^3, then 



— and this im- 



plies 



dxidxi 



U — 0. Therefore we have H^{x) = a^x\ + 03X1X3 and H2{x) 



4 I 22 

Qj^X-^ ~\- C4X1X2 



64X2- In the rest of this section we will not show all of the 



calculations. The next step to calculate the desired normal form is to transform 
Hi and H2 under the change of coordinates x = Pz, with help of a algebraic 
manipulator this is not problem, the interesting task is how to calculate the ho- 
mological equations ( 14 ) and ([15]) . There exist several ways on how to calculate 
these equations, in [3T], [53] and [TT] we can find a detailed discussions on this 
problem. Several authors use complex variables to perform the calculations of 
the normal form, however, it is worth noting that sometimes the normal forms 
under complex symplectic transformations do not determine the real normal 
forms [H] , therefore we use only real symplectic transformations to do the calcu- 
lations. In [3S] there is a elegant method to calculate the homological equations 
in the case non-semisimple of a linearized vector field and its application to the 
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R3BP at the equilibrium point L4, we are going to explain it briefly. 



Let us consider the change to symplectic polar coordinates {zi, Z2, Z3, Z4) — >■ 
(r, 9, R, e) given by 

Zi = r cos 9 
Z2 = r sin 6 

Z3 = R cos 9 — — sin 9 

r 

Z4 — R sin 9 + — cos 9 
r 

in these coordinates the term Hq becomes Hq = coQ + ^r^, by scaling time we 
can assume that = 1; therefore 

i/° = e + V (23) 

of course we must transform again the terms Hi, H2 under this change of 
coordinates. Let £,jjo{Wi) the linear operator defined by £^o(Wi) = {Wi,!!^}, 

i € N, explicitly ^HgiW,) = {W,, 9} + ^{W,,r^} = ^ -°rff := £5 + £iv 
with £s = and = ^''^^ The homological equation in every step looks 
like _ 

2ho{W,) + H', ^ H, 

where Hi groups all the terms of the previous steps Hj. Deprit ^11] proved that, 

after splitting Hi = Hi (r, — , i?, Q) + H^{i\ 9, R, 8), one can choose Hq = Hi 
and to solve 

ilHojW,)^Hl (24) 



with H'- = Hi — Hi in terms of operators we can write the equation ( 24 ) as 

in general the operator (id + £^^£jv) i^i) is not invertible, but in our particular 
case we can calculate Wi because all the terms of the right side of (24 1 are 
periodic in 9 and £,n is nilpotent. Because of the linearity of £,jjo we have 

{id + ils^£Ny^ --id- (S^s^S^n) + {^s^S,n)^ 

this series is a finite sum because in general 2IJ^{p) = if p is a homogeneous 
polynomial of degree m — 1. Therefore Wi is taken as 

W, = {id + Z-g^ilNy^ ^-s\H[) 

observe that we obtain a function totally periodic in 9. For this reason, we can 
repeat the process at any order. In particular after two steps we obtain 

G(6, y) = i?0(y) + eHliy) + '^Hl{y) 

The generating function can be known totally in every step without additional 
assumptions and it is entirely polynomial. Now we want to apply the method 
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explained above to calculate the Hamiltonian ( [T6| in the case of the restricted 
four-body problem. Let Hq be the Hamiltonian (|23|, as above, the first order 
homological equation to solve is 

If = Hf{r,-,R,e) + H^' {r,9,R,e) and we choose = Hf then the 
homological equation to solve is 



£ho{W,)^H'1 (25) 
with = — H'^ . In polar coordinates the term H'^ looks like 

ijj* — ^piir cos 9 + pu ^i? sin 6* + — cos ( 
C3 ( P22'f sin 6* + P23 ( i?cos 9 — — sin 6] ) ( pur cos 9 + pi4 ( i?sin 9 + — cos 9 



r J J \ \ r 

observe that the terms with 9 appear in all the terms of therefore — 

and = . If we write I — J H^d9 we can see that — 0, so the first 

order generating function is given by 

where the exponents stand by the composition of operators. For the second 
order term we have to calculate Hi = + {H^, Wi} + {H^, W2} and = 
Hi + {H^,Wi} but H^ = then H§ = Hi = H2 + {Hll,W2} with H2 = 
H2 + {Hi, Wi}, therefore the second order homological equation is 

ilHS{W2)=H^ 

with Hq — H2 and H2 — H^ + H'2, the second order homological equation 
is solved as before. Note that along this section we have performed all the 
calculations without using the numerical values of the coefficients of the matrix 
P and the polynomials Hi and H2. The expressions of Wi and Hq obtained with 
help of the algebraic manipulator are too long and we won't show them. Now if 
we substitute the required numerical coefficients of the matrix P: pu — —03928, 
Pi4 = -0.7631, P22 = -0.9680, P23 = -1.8807, P32 = 2.005, P33 = 1-3490, 
P4i = 0.8134, P44 — 0.5474 and of the polynomials ones 03 = —0.962, C3 = 1.370, 
a4 = -1.007, C4 3.150, 64 = -0.4686 we obtain 

H^ = hir-^ - h2R'^ + /i37-^i?^ - h^Q + h^R^Q - h^Q'^ - hT—:^+ (27) 

with hi = 2.19104, h2 = 1.41252, h^, = 16.2535, /14 = 8.35177, = 4.24874, 
/i6 = 0.00392, hj = 2.82504, = 4.24874, /ig = 1.41252. Therefore the 



Hamiltonian ( 16 1 after two steps looks like 



g = h'q + ^-hI 



11 



In the non-semisimple case the normal form is not unique but it depends on how 
we perform its construction, in [31] it can be found good discussion on this topic, 
however, a straightforward calculation shows that the superior order term (27) 
contains the terms of = {p\ +pl)[A{pl +pl) + B{qip2 - 92^1) + C{ql + gf)] 
and = c{qiP2 - q2Pi? + d(giP2 - 92Pi)(p? + pI) + e(pf + plf) with A, 
B, C, c, d, e constants, found in [55] and [53] in polar coordinates. This is 
a direct consequence of our construction of the normal form where we are not 
putting additional conditions on the generating function W and using only real 
coordinates. 



5 Analysis of the normal form 

We are interested in making a local analysis of the truncated system and of the 
stable and unstable manifolds. First, we will outline the construction of a versal 
deformation (|31|. [1]) of the normal form obtained in the previous section. It 
will useful recalling the definition of a versal deformation 

Definition 5.1. A deformation A{X) 0/ Aq e il/(C„) is called versal if any 
deformation B{pL) of Aq is equivalent to a deformation induced from A, i.e., 

Bi^,)^C{fi)Aim)C''M 
for some change of parameters : S — > A. 

This kind of deformation is useful to study the truncated system through 
resonance. We already said we do not have a explicit relation of the partial 
derivatives a and b with the mass parameter /i so we cannot give a explicit rela- 
tion between the parameters of the versal deformation with the mass parameter 
fj.. Therefore we just give a general relation of the parameters of the versal 
deformation with the quantities a and b. Consider the matrix B related to the 
linear vector field Hq, if we want to construct a versal deformation of the matrix 
B we must take the matrix family B^, := B + where the matrix Du lies on 
the orthogonal complement of the orbit of B [?], this is equivalent to ask that 
the matrix D^, satisfies = where _D* denotes the adjoin of Z?,^. The 

matrix = z/ici + ^262 with 



/O -1 



ei 



1 


V 



and 



62 



/ 



V 



1 






\ 



-1 
/ 



1 


/ 



satisfies is such that [D*,B] = 0. Consider the Hamiltonian matrix A of the 
linearised vector field of the R4BP 



A,. 



( ' 

-1 

a 

\ 



1 
1 

; 
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where a — i/aj^a;^ (0, /i) and Ux2X2{^t iA- We have found a hnear symplectic 
transformation P which brings A^-^ into its normal form B = P^^vl^^P, as we 
said the matrix 



V 





-1 




-(l + ^i) 



-1 



V2 




l + Vl 








(28) 



is a versal deformation of the matrix B. Consider now the family of matrices 
B^ = P^^A^P, B^ is a deformation of B (in /i) so, by the definition (5.1 1 there 
exists a change of parameters (/) : S — >■ A such that i?^ — C^D^C'^^ therefore 
and have the same characteristic polynomials, i.e., the polynomials Pb ^ — 
\^ + {2-a-b)\^ + {a+b+ah+l) and Ps„ = \^ +2[{l+vif + [{l+vif 
are equal. A straightforward calculation shows: 

V2 = 1/2 - l/4(a + 6) - l/2(a + b + ab+ if^ 

vi = [1/2 - l/4(a + 6) + l/2(a + b + ab+ l)i/2]i/2 _ i 



furthermore (1 + vi) 
expression 



1 — 1/2 (a + 6) — V2, the eigenvalues are given by the 
A' = -((1 + yif - V2) ± 2|1 + v^\^2 



note that the sign of V2 determines the behaviour of the eigenvalues i.e., whether 
they are purely imaginaries or with non-zero real part. In polar coordinates the 



versal deformation related to the matrix (|28|) looks like 

(1 + 



i?2 



92 



V2 2 



. The parameters v\ and V2 are not independent and is a first integral so we 
have that the versal deformation depends only of the detuning parameter V2. 
Now we are ready to make a study of the truncated system through resonance, 
the Hamiltonian under consideration will be 



2 



i?2 



62 



— 7' 

2 



e 



Hi 



Because is a constant we have obtained a system with one degree of freedom, 
we want to see which terms are the most important near the origin so we will 
use the scaling of (24j in polar form 

r ^tr (29) 



V 



which is symplectic with multiplier e The Hamiltonian becomes 
H^Q + e(\(R^ + ^\+Y + h,A + 0{^) 



(30) 
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If we ignore the higher order terms, we obtain the form of the truncated Hamil- 
tonian considered in |24j with parameters 6 = 1 and rj = hi > 0. So from the 
conclusions of that work, we have that the stable and unstable manifolds of the 
truncated system are globally connected for < and shrink to the equilibrium 
as — > 0~ , actually, when we consider small perturbations of the truncated sys- 
tem or equivalently, when we consider the higher order terms of the full system 
for e sufficiently small the stable and unstable manifold may not agree but they 
will still intersect. This last affirmation is a an application of the Poincare's 
argument that basically tells us that nondegenerate homoclinic points are sta- 
ble under small perturbations. In the references jH] and [20] the reader can 
find detailed discussions on these topics. Therefore, we have satisfied the re- 
quired conditions under which the "Blue sky catastrophe" termination occurs: 
a topologically transverse intersection of the stable and unstable manifolds of 
the equilibrium point L2. 



6 The unstable and stable manifolds 

In this section we perform a brief numerical study of the unstable manifold 
and the stable manifold W" to show the transverse intersections which were 
predicted in the former section. We do not attempt to make an exhaustive 
numerical study of such manifolds so we will restrict ourselves to show some 
transverse intersections for a few values of the parameter fi we were able to cal- 
culate numerically. We recall that from the linearization around the equilibrium 
point L2 was discussed in the section 3 we saw that for that point there exist 
a value of the mass parameter such that for all /z € {fib, 1/3] the eigenvalues 
given by ([9]) are Ai^2,3,4 = ± iuj where a and uj are real and strictly positive 
quantities. The general theory states that for this case the solution of the linear 
part of the equations ^ is given by 

4 

x{t)=Y.^^^^'' (31) 

i=l 



Ait 



the coefficients are not independent. In [30] it can be found that such coefficients 
satisfy the relation (Af — Q,xx)Ai = (2A,; -I- flxy)Bi but the equilibrium point is 

„ , A? - 



coUinear so = then 



2A, 



However, as we have seen along this paper, we can not know analytical relations 
between the eigenvalues and the mass parameter and therefore we do not have 
a explicit relation between the constants Ai ,Bi and the mass parameter /i. So 
for a numerical globalization of the invariant manifolds and we have 
decided to take a set of initial conditions totally equivalent to the classical initial 
conditions shown in T4 . Our vector of initial conditions is given by 

= L2 + e(cos 6v + sin 6w) (32) 
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where L2 — {xl2,070,0), e > 0, 6 & [0, 2n] and v and w are in the unstable (sta- 
ble) space i?" {E''). As we have already mentioned, the equations of motion have 
the property that if a; = x{t), y — y{i) is a solution, then x = x{—t), y = —y{—t) 
is also a solution so the unstable and stable manifolds of the equilibrium point 
are symmetric with respect to the x— axis and therefore is enough to calculate 
the manifold VF". In the R3BP and more recently in the restricted four body 
problem (where the primaries are in a coUinear configuration) several authors 
[Hj, [5] have studied the "Blue Sky Catastrophe" termination and the asymp- 
totic orbits by analyzing the transverse intersections of the invariant manifolds. 
In these cases the equilibrium points, where the BSC termination ocurrs, are off 
the a;— axis but in our case the equilibrium point L2 is on this axis, so we have 
to be careful when developing the numerical globalization. As in the R3BP we 
use the {x, x) plane to represent the transverse cuts of the invariants manifolds, 
this can be done because we are using the surface section y = and because 
the velocity y can be put as function of x and x from the relation Using ^ 
and ( 32 ) we compute and show in the figure ([7| some transverse intersections 
of and on the symmetry axis. From the linear approximation we have 
that the fiow is of spiral type in a small neighbourhood of the equilibrium point 
L2 then the orbits on the unstable manifold have several intersections (cuts) 
with the surface section y = 0. Therefore in the "first" cuts we do not expect 
having a transverse intersection i = of in the plane {x,x), so we have 
to consider the cuts when the orbits leave an appropriate neighbourhood that 
depends of the parameter e. In the figure Q we can see that after four cuts 



m3 




"^2 





Figure 4: The fourth cut of the unstable manifold of L2 in the (x, y) plane for 
= 0.019. 

the orbits do not have orthogonal intersections (i = 0) with the symmetry axis 
and therefore the unstable manifold does not intersect the x— axis in the (x, x) 
plane, when we continue the integration we observe that the infiuence of the 
primaries 1122 and 7713 appears on the orbits, if we use a regularization process 
this is not problem, the interesting phenomena is that the continuity of the 
observed in the first cuts is lost. In the figures ^ and ^ we show the unstable 
manifold for /i = 0.2, in order to avoid congestion of the picture we present the 
manifold for e [0, tt] and 9 e [tt, 27r]. 

However we can find transverse intersections of on the symmetry axis. 
For the value of the parameter /i = 0.019 we found the transversal intersection 
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Figure 5: First row: the fourth cut of VF" for /i = 0.2 and 9 G [0, tt]. Second 
row: magnifications of the intersections with the x-axis 
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^'^^^^ forbidden area / 
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Figure 6: The fourth cut of W for ^ = 0.2 and 9 E [n, 2tt] 



P2 ~ 1.925 corresponding to the honiochnic orbit that is the termination of the 
family / see figure If the value of /i is increased, more distant transverse 
intersections of the three primaries appear, see figure ([t]). 
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Figure 8: The blue sky catastrophe termination for the family /. 
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